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Abstract 

We consider the e-regime of QCD in 3 dimensions. It is shown that the leading term of the 
effective partition function satisfies a set of Toda lattice equations, recursive in the number of 
flavors. Taking the replica limit of these Toda equations allows us to derive the microscopic spectral 
correlation functions for the QCD Dirac operator in 3 dimensions. For an even number of flavors 
we reproduce known results derived using other techniques. In the case of an odd number of flavors 
the theory has a severe sign problem, and we obtain previously unknown microscopic spectral 
correlation functions. 
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1 Introduction 



Recently there has been substantial progress in understanding the so-called replica method in the 
context of the effective low-energy field theory for QCD (QCD4) [1-7]. While perturbative (series) 
expansions readily yield correct results with this method [1], it has proved difficult to obtain exact 
nonperturbative results [8, 9, 10, 11]. A crucial ingredient in formulating a systematic replica method 
which is reliable also for nonperturbative calculations has been the use of exact relations of the Painleve 
[4] or Toda [5, 6, 7] kinds. These relations reflect the fact [12] that the leading-order effective partition 
function in the so-called e-regime [12] 1 of QCD is what is known as a r-function of an underlying 
integrable KP hierarchy [14]. This integrable structure derives back to the spontaneous breaking of 
chiral symmetry according to the pattern SUl(Aj)xSUr(Aj) — ► SU(Aj), where Nf is the number of 
light quark flavors. As a consequence of this spontaneous breaking of chiral symmetry and because 
of the mass gap in QCD the low-energy dynamics of QCD is governed by the Goldstone modes. 
In the standard counting scheme the low energy theory describing the dynamics of the Goldstone 
modes is known as chiral perturbation theory, and in this framework the ordinary replica method 
applies without subtleties [1]. If one instead [12] considers a counting scheme (the e-regime) where the 
Compton wavelength of the pions is much larger than the linear dimension of the volume, l/m n S> L, 
then the partition function at leading order reduces to a static integral over the Goldstone manifold. 
Changing Nf in this integral builds up a series of partition functions which are connected through a 
Toda lattice equation. This connection between partition functions with different Nf shows how to 
take the replica limit Nf — ► 0. The purpose of this paper is to establish Toda lattice equations for 
QCD in 3 dimensions (QCD3) and to show that also these allow us to obtain exact non-perturbative 
results using the replica method. 

QCD3 cannot undergo spontaneous breaking of chiral symmetry in the usual sense. Moreover, also 
the notion of gauge field topology is different, a remnant manifestation being the possibility of adding 
a Chern-Simons term to the action. Nevertheless, a three-dimensional analogue of chiral symmetry 
breaking (really flavor symmetry breaking, re-interpreted) is possible [15]. The basic idea is that 
the chiral components of the four-dimensional spinors in three dimensions can be mimicked by two 
different fermion fields, with masses of equal magnitude but opposite signs. Such mass terms of 
different signs are possible in an odd number of dimensions, a consequence of the fact the two sets 
of 7-matrices {7i} and {— 7i} form two inequivalent irreducible representations of the Clifford algebra 
in an odd number of space-time dimensions. The two 2-spinors corresponding to opposite signs of 
the mass terms may be grouped into one four-spinor of which the top components play the role of 
the left-handed field, while the lower components represent the analogue of the right-handed field 
in four-dimensional language. Moreover, the associated "chiral" symmetry can break spontaneously, 
and is expected to break spontaneously when the number of fermions species is small. The symmetry 
breaking pattern is believed to be that of [15, 16] 

U(2N f ) -+ U(N f ) x U(Nf). (1.1) 

There are numerical simulations of QCD3 on rather small lattices that support this conclusion in the 
extreme (quenched) case of no dynamical quarks [17]. An odd number of flavors is obviously a more 
difficult issue. For 2Nj + 1 flavors of which the 2Nf masses are grouped into pairs of equal magnitude 
but opposite signs it has been argued in ref. [16] that the spontaneous symmetry breaking pattern 
(1.1) is replaced by 

U(2N f + 1) U (Nf + 1) x U(Nf) . (1.2) 
1 See e.g. ref. [13] for a recent review with an emphasis on the present context. 
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This case is rather tricky to treat from the effective field theory point of view. An interesting starting 
point is Random Matrix Theory, which suggests that there are two possible effective field theories, 
depending on the number of Dirac operator eigenvalues [16]. We will look into this problem in closer 
detail below. We note that as one further generalization one could also consider the case of 2Nf flavors 
paired with opposite signs, and n unpaired flavors. 

The Random Matrix Theory representation of the effective partition function for QCD3 in the e- 
regime [16] has lead to a number of intriguing conjectures. In particular, spectral properties of the 
pertinent 3-dimensional Dirac operator have been argued to be derivable from the joint eigenvalue 
probability distributions of such a Random Matrix Theory [16, 18], and based on this conjecture 
universal spectral n-point correlation functions have been found for both an even [19, 20, 21] and 
an odd [22] number of flavors. First steps toward deriving these results directly from the effective 
field theory partition function were taken in refs. [23, 3]. In particular, by an impressive series of 
supermanifold integrations some of the spectral correlation functions were derived from the so-called 
"supersymmetric" formulation in ref. [23]. In [3] intriguing relations between the QCD4 and the QCD3 
partition functions in the e-regime were derived. Applying the replica method to these relations a set 
of identities between the spectral correlation followed. These identities, derived by formally applying 
the replica limit Nf — > 0, presupposed a meaningful operational way to perform this limit, a non-trivial 
issue in view of the known subtleties involved [8, 11]. Here, using the replica limit of the Toda lattice 
equation, we show that these results, and others which are new, can be given a precise meaning. 

The new results are for an odd number of flavors. This theory is rather special since after integrating 
out the fermions the weight inside the partition function is not positive definite. This manifest itself 
in the eigenvalue density in a number of ways. For example, it is tempting to conclude that eigenvalue 
density in the theory with and odd number of flavors and an odd number of Dirac operator eigenvalues 
[16] is an odd function in the chiral limit. However, as we will explain below this is not the case. 
Although the non-positive weight in QCD3 with an odd number of flavors comes about quite differently 
from the sign problem in QCD4 with a non-zero baryon chemical potential, the lessons learned from 
QCD3 may be worth keeping in mind when tackling the sign problem in QCD4. 

The presentation of this paper is as follows. First we derive the Toda lattice equations for QCD3 
in the e-regime. Then we discuss the extension to bosonic flavors. Given the Toda lattice equations 
we take the replica limit in section 3 and show how to obtain spectral correlation functions. As in 
QCD4, the Toda lattice equations in QCD3 also hold for partition functions with both fermions and 
bosons. This allows us to establish the replica limit of the Toda lattice equation in QCD3 on the same 
footing as the supersymmetric method [24]. We show how in section 4. Finally, in section 5 we show 
how general Toda lattice equations follow from one single Consistency Condition of Random Matrix 
Theory. We sumarize our results in section 6. 

2 The Toda Lattice: From QCD 4 to QCD 3 

Our starting point is a set of Toda lattice equations which have been derived [25, 26, 14, 6, 7] for 
the leading term of the effective QCD4 partition function in the e-regime. Using these and known 
connections between QCD4 and QCD3 in the e-regime we establish the Toda lattice equations for 
QCD 3 . 

We first give our conventions and definitions, which follow ref. [3]. The e-regime is in this (2+1)- 
dimensional context defined by the large- volume limit where nevertheless V <C \jm\ and m n generi- 
cally denotes the pseudo-Goldstone masses of spontaneous breaking of chiral symmetry. In this "ex- 
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treme" chiral limit where m = niiVE is kept fixed the partition function reduces to a static integral 
over the Goldstone manifold. The infinite-volume condensate is defined by 



£ = lim lim , (2-1) 

m~*0 V— >oo 



and rrii are the quark masses. 



The leading contribution to the QCD3 partition function in the e-regime is for an even number 2Nf 
of fermions grouped in pairs with masses of equal magnitude but opposite signs given by [16] 

z ? Nf) ({M) = / dl/exp[VERe Tr{MUT b tf)] , (2.2) 

JU(2Nf) 

with Ai=diag(mi, ... ,miy f ,— mi, ... ,—miy f ), T5 = diag(lAr / ,-^-N f ) and dU always indicating the 
Haar measure. This integral can be performed explicitly with the help of the Itzykson-Zuber formula 
[27], resulting in [20] 

det / A({^}) A({-^}) \ 

*T'\M) = i-iY«iw ^ A({ "^) A({M) ] ^ 

where, with our normalization conventions, the Nj x#j matrix A({^ij}) is defined by 

MMh = (H) 1 ' 1 ^ > j,l = l,...,N f , (2.4) 
and the Vandermonde determinant is given by 

2N f 

MM) = HbH-Hj) ■ (2.5) 

i>j 

As mentioned in the Introduction, the case of an odd number of flavors is a bit unusual because the sym- 
metry of the action relies on a paring of masses. Nevertheless, as in the even-flavored theory, the effec- 
tive partition function has a unique leading-order term in the e-scheme. With .M=diag({//j}, fi, {—fJ-i}) 
and T$ = (ljVf+i, — ljV/) the partition function becomes [16], 

zS N,+1 \fi,{fi i }) = [ (f[/cosh[V£Re Tv(MUF 5 U r )} , (2.6) 

JU(2Nf+\) 

for an even number of Dirac operator eigenvalues, and 

Z^ Nf+1 \n,{ni}) = f dUsmh[VXRe Tr(MUT 5 tf)] , (2.7) 

JU(2Nf+l) 

for an odd number of Dirac operator eigenvalues. Both of these integrals can again be evaluated 
explicitly, with the conventions of [3], 

Z<f' +1) (MM) = (-l)^^ +3 )/ 2 ^y-[detD( M ,{^})±(-l)^detD(-^{- W })] , (2.8) 

where the (2Nf + 1) x (2Nf + 1) matrix D is defined as 

2^detD(MM) ee detf A (MMV,+i*at /+ i A(- M {-MK+ixiv, \ ( j 

\ A({-m}) NfXNf+1 A({fn}) NfXNf 
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and A is denned in eq. (2.4). 

For QCD4 in the corresponding e-regime the leading contribution to the partition function is 2 

Zv Nf \{»i}) = I dU (detU) u exp[VT,BjeTr(MU)] (2.10) 

Ju(N f ) 

in a sector of topological charge v [12, 28]. The quark mass matrix in 4 dimensions is A4=diag(mi, . . . ,rriN f ). 
Also this integral can be explicitly evaluated for integer v [29]: 

where the matrix B in eq. (2.11) is given by 

B({ft})ii = ^ 1 4 i_1) (^) , j,l = l,...,N f , (2.12) 
with I® the /'th derivative modified Bessel function I u , and the denominator is given by the Vander- 



monde determinant of, in this case, squared quark masses, 

{A - m?) = 



A(K 2 }) = n^ 2 -^ 2 ) = det {fir' 1 ■ (2-13) 

i>j 



We use on purpose the same notation ^ in both QCD3 and QCD4, since relations between these two 
widely different theories exist when these dimensionless parameters are identified in the two theories. 
The QCD4 partition function for non-integer u is defined by analytical continuation in the index v of 
the Bessel functions [3]. We stress at this point that the explicit formulas (2.3), (2.8) and (2.11) of 
course are only valid for Nf taking positive integer values. Below we shall return to the issue of how 
to deal with cases where Nf is zero, or even negative. 

2.1 Two theorems 

Our main results will be based on the following two relations between the QCD3 and QCD4 partition 
functions. These relations were proven in ref. [3] and hold when the normalization conventions are as 
stated above. 

Theorem I 

4 2Nf \M) = ^42 /2 ({M) 4!i /2 ({^}) • (2.i4) 



Theorem II 

= ^'V^*2$(MM) 4S /2 ({M) • (2.15) 

We have written these theorems explicitly for real and positive masses {/Xj}. When taking discon- 
tinuities in the complex plane one should always recall that only the absolute value enters here. In 
what follows we will for notational simplicity normally not explicitly display the associated factors 
that equal unity for real and positive values (usually terms of the form x / \fx^) . 

2 We use a notation in which we do not explicitly indicate that this is the 4-dimensional partition function; the partition 
function is seen to be the one for QCD4 by the labeling according to topological charge v. There should be no source of 
confusion as we will never explicitly consider the case v = 3. 
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2.2 Toda lattice equations for QCD 3 



The two theorems (2.14) and (2.15) can nicely be combined to derive Toda lattice equations for QCD3 
on the basis of the known Toda equations for QCD4. As a first illustration of this procedure, consider 
the case of Nf degenerate flavors of mass x in QCD4. In this case the leading term of the QCD4 
partition function satisfies the differential equation [25] 



(xd x ) 2 lnZ { u Nf \x) = 2N f x 



2 %v {x^^v (^0 



(2.16) 



We can turn this into a Toda equation for QCD3 by use of Theorem I and II. By considering the case 
of an even number of flavors 2Nf with degenerate masses x we get 



{xd x ) 2 \nZf N! \ X )=AN f x^ 



Z^ +1 \x)zf N ^\x) + Z^ +1 \x)Z^~ 1 \x) 



(2.17) 



This is the QCD3 Toda lattice equation analogous to the QCD4 equation (2.16). We note the quite 
generic feature of the QCD3 equation having two contributions on the right hand side, originating 
from having two different partition functions with an odd number of flavors. Of course, we might have 
derived the Toda equation (2.17) directly from the explicit formulas (2.3) and (2.8). But the present 
shortcut through use of Theorems I and II is obviously a much simpler route. 

We next consider the QCD4 Toda lattice equations with Nf degenerate flavors of mass x and n 
degenerate flavors of mass y. These are [25] 



and [6] 



xd x (xd x +yd y )lnZ^ n \x,y) = 2N f x 2 ^ r - 2 , 

zl Nf ' n \x,y) 

zi N - n \x,y)} 2 



(2.18) 



(xd x yd y )lnZl Nf ' n \x,y) = 4N f nx 2 y 2 



(2.19) 



The corresponding Toda lattice equations for the QCD3 partition functions are readily obtained from 
the theorems. To keep it simple, we will start by looking at 2Nf flavors with degenerate masses ±x 
and 2n flavors with degenerate masses ±y (both with the usual pairing). After re-expressing the Toda 
equations in terms of QCD3 partition functions we get 



xd x {xd x + yd y )\nZl f ' '{x,y) 

Once again the right hand side has two contributions. The generalization to the case where the 2n 
masses are paired but non-degenerate is obtained simply by replacing y with {yi} and (xd x + yd y ) by 
(xd x + J2?=i Vidyi) ■ Following the same route as described above, we find that the QCD3 equivalent 
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of (2.19) for an even number of flavors is 



xd x yd y log zf Nf ,2n) (x, y) 



(2.21) 



|V , , ' zg^ +2 ' 2 " +1) (*, y)z£ N < ^ (x, y) Z £ N ' + ^ (x, y)Z^ ~ 2 ^ ( X ,y) 



z<r ,+1,an W)* 



(2N f -l,2n) 



^ +1 ' an) (x,y)^'- 1 ' aB) (x,y) 



Despite appearences, the right hand side is symmetric under interchange of x and y. 

We now turn to the partition functions with an odd number of flavors. Starting at the simplest case 
with 2Nf flavors of paired masses ±x and one additional flavor of mass x we get from (2.16): 



(xd x ) 2 \nZ ( p +1) (x) 



(2.22) 



= 4x 



(N f + iy. 



zf ± N ^\x) Z f^ +l \ 



V 



z. 



(2N f +2) 



(X) 



+ N f - 



7 {2N f +l) 
J 3* ' 



x)Z. 



(2N f -l) 

3± 



z. 



(2N f ) 



(X) 



For the QCD3 version of (2.18) where we have 2Nf paired and degenerate masses ±x and 2n + 1 
flavors of mass y we find 



xd x (xd x + yd y ) In zf± f,2n+r> ( x , y) 



= 4N f x 



7 (2N f +l, 
■-3± 



2n+2 \x,y)Z^ 2n+2 \x,y) , Z^ 2n \x,y)Z^^x,y) 



+ 



4 2Nf > 2n \*,y) 



[zf N - 2n+2 \x,y)] 

As the final Toda lattice equation for QCD3 we also give the analogue of (2.22), 

xd x yd y log zf± !,2n+1 \x, y) 

„3„, I 1 i\^3± V^'i/V^T 



4 2 ± ^ +2 ' 2 " +3) (x, y )4f^ 2 ' 2n+1) (x, y ) 



= 47V / x a y[(n + l) g 2 ± JV/+12n+2) 

^ /W i) (g>y) ^ 



+n- 



(*,y)^'- W2 W) 
x,y) 



3± 

(2JV/-2,2n-l). 



^ +1 ^(x,y)^'- 1 ^(x,y) 



(2.23) 



(2.24) 



valid for an odd number of flavors. 



2.3 Extension to bosonic flavors 

Before proceeding to use the replica method we must ascertain that the Theorems I-II are consistent 
with what has been understood about the QCD 4 partition functions for both zero and a negative 
number of flavors. We should first clarify what is meant by this. For zero flavors the partition functions 
are "quenched", the fermion determinants are entirely absent, and there can be no mass dependence. 
It is clearly convenient to choose the normalization so that in this case the partition function simply 
equals unity. A negative number of fermionic flavors is defined by raising the determinant to the 
corresponding negative number. This is equivalent to a partition function of the same number of 
complex fields with bosonic statistics (and thus violating the spin-statistics theorem, but these bosons 
are never considered as external physical states) . On the QCD4 side the generalization of the partition 
function formula (2.11) was given in ref. [30, 5, 31] for an arbitrary number Nf and Nf, of fermionic 
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and bosonic species, respectively. For positive Nf and Nj, the explicit expression reads, in a hopefully- 
obvious notation, 

zi NANb \{x f }\{y b })= ^ lj ^-^=^f^ , (2 . 25) 

where = Xj for i = l,...,N f , z Nf+i = yi for i = l,...,N b , J u +j-\{zi) = I v+j -i(xi) for i = 
1, . . . , Nf, and J v +j-i{zN f +i) = (— for z = 1, . . . , TV;,. The generalization to negative 

integers is done through identifications such as zi Nf ^ Nb ^ = £^ Nf+Nb ^ and so on. Thus, while the 
Nf = 1 QCD4 partition function is given by Z^\x) = I w (x), the Nf = — 1 partition function reads 

= K v {x), where K u (x) is the modified Bessel function. We therefore know how to define the 
right hand side of the two theorems for a negative number of flavors. The question is whether we can 
define the left hand side for zero or a negative number of flavors in QCD3 in this way. Setting Nf = 
in eq. (2.14) gives Z^ = 1 from Theorem I, i.e, with the correct normalization. Having gained some 
faith in this procedure, we can infer the QCD3 partition function for Nf = — 1 from Theorem II: 

Z^\x) = {^K ±l/2 {x) = e~' • (2.26) 

This agrees exactly with what we should expect on general principle: For Nf = — 1 (one boson) there 
is no spontaneous symmetry breaking at all, and the leading term of the partition function is simply, 
by a generalization of the argument by Leutwyler and Smilga in the 4-dimensional case [28], the 
exponentiation of the leading term free energy F = mVT, = fj,. It thus appears that we can continue 
with the indicated identifications, and Theorem I will then give us the partition function for Nf = —2: 

Z { f 2 \x) = -K_ 1/2 (x)K 1/2 (x) . (2.27) 
This procedure obviously continues for higher (negative) values of Nf. 

We can formalize the above arguments by making the following conjectures regarding supersymmetric 
generalizations of Theorems I and II: 

zi 2N ^\{x i} m) = ^/-^^^({Xi}!^})^'^^}!^}) (2.28) 

zii N ^\ { x,}\{y l} ,y) = ^-^yi^^^xOK^},,)^^^}!!^}) .(2-30) 

In all cases we have been able to check these identities are consistent with what we know from other 
sources. We note in particular that if they are correct we have circumvented the quite difficult task 
of evaluating all the so-called Efetov-Wegner terms in the supersymmetric version of the effective 
Lagrangian [23]. It appears that the most direct way to prove these identities may be through an 
explicit evaluation of the related supersymmetric Random Matrix Theory integral, as in ref. [31]. 



3 Replica limit of the QCD3 Toda lattice equations 

As an introduction let us consider QCD4 with Nf quarks with masses {/x/} and n fermions of mass x. 
Using the replica method the partially quenched chiral condensate, the resolvent, is then defined as 

Gi Nt \x,{pL f }) = liml^ln^' n) (W},x) . (3.31) 
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The spectral density of the Dirac operator in QCD4 with Nf flavors is given by 

pl N '\\,{» f }) = (j25(\ k -\)) , (3.32) 



where (. . .)N f is the vacuum expectation value in QCD4 with Nf flavors. The delta functions can be 
expressed as the discontinuity of the resolvent across the imaginary axis 



pl Nf \\,{»f}) = ^DiscGf^z, {p f }) 



x=i\ 



= ^-hm[Gi Nf \iX + e)-G { l Nj) (iX-e)} . (3.33) 

The resolvent and the density in QCD3 with an even number of flavors 2Nf follow analogously [3]. 
In the replica formulation the resolvent is defined from the partition functions with 2Nf flavors with 
paired masses {±/U/} as well as 2n replica flavors with paired masses ±x: 

G^M^f}) = limi-A lnZ (-/.-) ({/X/} , x) . (3.34) 

As for the partition functions we use a similar notation as in QCD4, the two are easily separated by 
the explicit topological index v. 

The eigenvalue density is defined as in (3.32) except that the average is in QCD3 with 2Nf flavors. 
However, because the replicated flavors have paired masses the discontinuity over the imaginary axis 
in the complex quark mass plane of the resolvent becomes 



i-Disc Gf^*, W}) _ = ^te(S »rT^T3 + 



1 2e 2e \ 

2T72 / 

2N 



f 



x=i\ 4vr e-o \ Y (A fc + A) 2 + e 2 (X k - A) 2 + e 2 

= \{p? Nf \-\M)+P? Nf) ^{^f})) ■ (3-35) 

The limit n — > in the defining equations for the resolvents (3.31) and (3.34) must obviously be taken 
with care, since the partition functions entering the right hand side are only known for integer values 
of n. For more than two decades it was widely believed [8] that one could at best obtain small or large 
argument expansions of the true result using the replica method. Recent developments, see e.g. [9, 10], 
attempted to go beyond this, but also that approach was met with some criticism [11]. However, with 
[4, 5, 6, 7] this situation has drastically changed. One tool to perform the replica limit correctly is 
the Toda lattice equations for the leading-order QCD3 partition functions. How it works in detail is 
perhaps best explained by working out a couple of examples. 

Let us start with the simplest case in QCD3 namely the fully quenched spectral density of the Dirac 
operator. In order to obtain this we first determine the fully quenched resolvent from the Toda lattice 
equation (2.17) 

JLim ^( xd x) 2 log Zf n) (x) = xd x xGf> (x) 

= 2x (Z&> (x)Z { -- 1) (x) + (x)Z^ (x)) , (3.36) 

where we have used that Z^ Nf °' ) = 1. The partition functions can conveniently be found from (2.15) 
which then gives 

d x xGf\x) = xK 1/2 (x)(I_ 1/2 (x) + I 1/2 (x)) = 1 . (3.37) 
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The solution with boundary condition Gf^ (x) finite in the limit x —> follows readily, 

Gf\x) = l. (3.38) 
Taking the discontinuity of the resolvent we get 

^(pf(-A) + /4 0) (A)) = 4 0) (A) = l (3-39) 

where we have used that the even-flavored quenched spectral density is an even function in A. Eq. 
(3.39) is just the result derived from Random Matrix Theory [16], where, due to the absence of the 
determinant, it simply corresponds to the bulk microscopic UE density, which is flat. 

Next, we work out the eigenvalue density in QCD3 with an even number of flavors. Using the same 
technical steps as in [5] the resolvent in QCD3 with an even number of flavors follows from the Toda 
lattice equation (2.20). Consider the easiest example Nf = 1 which gives 



^0 2n 



(3.40) 

where once again the partition functions on the right hand side can be found from (2.14) and (2.15). 
This can be integrated to give a resolvent 

(2). . x 2 - y 2 - x + ye~ x (coth(y) sinh(x) + cosh(x) tanh(y)) , . 

and from this the density in QCD3 with two paired masses follows: 

( 2 ) ( 2 ), , 1 j/(cos 2 A tanh y + sin 2 A coth y) 

p\'{\y) = p\\-\y) = n(\ 2 + y 2 ) • ( } 

This is in complete agreement with the result found in Random Matrix Theory [20]. 

Now we move on to the odd sector. We look at (2.24) with one real flavor and 2n replicas for both an 
even and odd number of eigenvalues. Taking the replica limit of this Toda lattice equation gives us 

l^xdxixdx + ydy^nZ^^y) = 2x \ ^^k^^l +Z W(x)z£ L \z) | . (3.43) 



4\y) 



We start by looking at the case of an even number of eigenvalues. In this case we get the following 
resolvent 

(i) _ x 2 - y 2 + e~ x (y cosh(x) tanh(y) - x sinh(x)) 

^3+ { x > y) — x 2 -y 2 ' \ 6 -^) 

which then gives the symmetric part of the eigenvalue density in QCD3 with one flavor and an even 
number of eigenvalues 

\ (p« <-*. »)+/#(*. »)) = \ - AcosAsi :( A ti7) Ata "" i ' ■ < 3 - 45) 

In the massless case, y = 0, this agrees with the prediction from Random Matrix Theory [22] and 
thus gives indirect confirmation of the conjectured form of the supersymmetric partition functions in 
QCD 3 . 
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We next consider the case of an odd number of Dirac operator eigenvalues, for which microscopic 
spectral correlation functions have not been evaluated previously due to the unusual behavior of the 
partition function. This behavior does not obstruct our replica approach, in fact the evaluation of 
the spectral density is completely analogous to that in the even sector. The replica limit of the Toda 
lattice equation in this case gives the resolvent 

G (i) = x 2 - y 2 + e- x {y sinh(x) coth(y) - x cosh(x)) 

^ x 2 — y 2 

which then leads to the even part of the density in QCD3 with one flavor and an odd number of 
eigenvalues 

1 / (i) (i) \ 1 A cos A sin A - y sin 2 A coth y 
-^(-A, y )+^(A,y)j = - + 7T{X 2 + y2 ) • ( 3 - 47 ) 

This new result is plotted in figures 1 and 2, the plots showing respectively the massless case (y = 0) 
and a massive case (y = 10). Note in particular that in the massless case the sum does not vanish: 
the eigenvalue density in the massless case is not an odd function of A. Let us clarify this point. It is 
tempting to assume that the spectral density of the Dirac operator eigenvalues will be odd in A. Our 
explicit calculation above shows that it is not the case, but one can gain a better understanding of 
this phenomenon if one is willing to use the Random Matrix Theory representation for the eigenvalue 
density [16]. To this end, let Z denote the partition function in the Random Matrix Theory (only in 
the microscopic limit does this theory agree with the static field integral [16]), 

„ N N N N f 

4 N ii\{m f }) = \\d\ k n iA / -A fc i 2 ne- w(A ^ ) n(^+^/) ^ 
j k=i k<i=i k=i f=i 

for a suitable potential V(X 2 ) [19]. We will use this representation to highlight the special properties 
of the density in the theory with both Nf and N odd. It is sufficient to consider Nf = 1. First note 
that the partition function in this theory goes to zero linearly with m. This in itself is not alarming, 
and for example also the effective QCD4 partition function (2.11) vanishes like m u in the limit m — > 0. 
However in that case the analogous Random Matrix Theory includes an explicit factor m v in front 
of the eigenvalue representation, and this is easily canceled out when calculating the spectral density 
from the matrix integral. The small-mass behavior of the integral (3.48) is more complicated, as 
the mass m mixes with the eigenvalues. To be precise, we simply define the spectral density by the 
expectation value (3.32). From eq. (3.48) it then follows that the eigenvalue density is given by 



, (Ai + im)e-wW J UL2 dX k UH<i=i \>H - X k \ 2 IltU e~ NV ^(X k + im) 
p 3>N {M,m) = (1) . (3.49) 

Z K '{m) 



For odd N and Ai 7^ the eigenvalue density diverges like 1/m for m — > 0. However, for Ai = the 
m — > limit is finite. Using the eigenvalue representation above it is also easy to show that 

P^ N (-X, m) = p ( ^ N (X, -m) (3.50) 

and that the sum 

p { ^ N (X,m) +p { ^ N (-\,m) 

appearing in (3.47), is finite for all values of A even in the limit m — > 0. The divergence as m — > 
thus resides entirely in the odd part of the eigenvalue density, and this odd part is not probed by the 
discontinuity of the resolvent. 
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Our results in the odd-Nf theory with an odd number of Dirac operator eigenvalues therefore seem 
to be well understood also from the Random Matrix Theory representation. As one additional non- 
trivial check on our results, we note that the eigenvalue density (3.47) satisfies the correct decoupling 
condition of approaching the fully quenched spectral density when the mass y goes to infinity. 



0.5r \ (pW{-\,y)+pV>{\y)) 
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Figure 1: Plot of the symmetric part of the microscopic spectral density in QCD3 for one massless 
quark and an odd number of Dirac eigenvalues. The spectral density takes the value 1/tt at A = 0. 
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Figure 2: Plot of the symmetric part of the spectral density for one relatively heavy quark, y = 10, and 
an odd number of Dirac eigenvalues. The amplitude of the oscillations is smaller than in the massless 
case. As the mass, y, goes to infinity the amplitude decrease as 1/y, and the density converges toward 
the quenched spectral density, 1/tt, as expected. 



So far we have given several examples of how the eigenvalue density can be derived from the replica 
limit of the Toda lattice equation. One can also derive the chiral susceptibility in this fashion. This is 
in fact much simpler since the Toda lattice equation already has a double derivative, and no integration 
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is therefore required. 

To determine the chiral susceptibility 

Xz\x,y) = lim -^—^-^-\nZf m,2n) (x,y) , (3.51) 
n->-o,m->o 4nm ox ay 

we take the replica limit of the Toda lattice equation, (2.22), and find 

Using Theorem II and (2.30) we find 

xf(xy) = £ T!*l). (3 . 53) 

The double discontinuity of this function gives the double symmetric combination of the two point 
function [3] 

p(°)(A 1 ,A 2 ) + 4 0) (-A 1 ,A 2 ) + pf (A!,-A 2 )+pf (-A 1 ,-A 2 ) = -i- I Disc X f (x,y) x x (3.54) 



x=i\i,y=i\2 



We note that all of these cases, and any other we can consider, give the detailed support for the general 
relations that formally can be derived on the basis of Theorems I and II and the replica definitions 
(3.31) and (3.34) alone [3]. These are, for the one-point functions, 



for an even number of flavors, and 

p^ +1) (A;{^},/i)+4 2 + ^ +1) (-A;{/i,}^) = pS 1 /2(A;{^},M)+pS 1/2 (A;{/i,}) , (3-56) 
for an odd number of flavors in the "+"-sector. 3 

In ref. [3] the case of an odd number of flavors in the "-"-sector was not considered, but we can easily 
derive the analogous relation by use of Theorem II. We find 

pf Nf+1) (A; M, M ) +pf N f +1 \-\; {^}, M ) = p { u N J$ 2 (\-,M,») +pl N jl 1/2 (X;M) , (3.57) 

for which our (3.47) indeed is a special case. 

We also remind that similar general relations can be worked out for higher /c-point spectral correlation 
functions [3], and that Toda lattice equations can give the justification for the formal replica manipu- 
lations that lead to these relations. Because of the pairing of masses with different signs these relations 
do not just involve the /c-point spectral correlators themselves, but combinations with different signs 
(see ref. [3] for details). We have already seen examples of this phenomenon for the spectral 1-point 
functions derived above. Also the chiral susceptibility (3.53) gives the symmetric combination of the 
spectral 2-point function. To isolate the spectral correlation functions directly one needs to combine 
the replica method with more than just the discontinuity formula of the mass-paired resolvent and 
higher-order versions thereof. 



1 



3 Note that ref. [3] incorrectly assumed that p^ f+1 \\; {fJ,i}, n) is an even function of A also in the case of fi ^ 0. 
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4 Graded Toda lattice equation 



While the Toda lattice equation allows us to take the replica limit in a well defined and correct way 
it does not fully explain why the analytic continuation in the number of flavors works. However, it is 
not difficult to understand why the replica limit of the Toda lattice equation and the supersymmetric 
[24] method give the same answer. To make this transparent [7] we consider supersymmetric versions 
of the Toda lattice equations and quench these as in the supersymmetric method [24]. Referring to 
the graded symmetry such Toda lattice equations have been called graded Toda lattice equations [7]. 

Two graded Toda lattice equations have been found for QCD4 [7]. Here we will take a closer look 
at one of these namely the supersymmetric variant of (2.18) with Nf degenerate fermionic flavors of 
paired masses ±x and Nj, degenerate bosonic flavors of paired masses ±y 

(N I AM Z {Nf + llNb) (x\v)Z (Nf - llNb) (x\v) 

xd x (xd x + yd y )lnZl N ^\x\y) = 2N f x^ r ™f - 2 ™ V) • (4.58) 

The generalized versions of Theorems I and II can be used to translate this equation in to a graded 
Toda lattice equation for QCD 3 

xd x (xd x + yd y )lnZi 2Nfl2Nb \x\y) 
m Z^ Nb \x\y^ 

fX r Zr^ Nb \x\y) 2 

It suffices to look at the special case of Nf = 1 and Nf, = 1, and quench it supersymmetrically. We 
start by focusing on the lhs 

lim \xd x (xd x + yd y )\nzf ],2 \x, -x\y, -y) = xd x xG$\x), (4.60) 

where we have used the supersymmetric definition of the fully quenched resolvent 

Gf\x) = hmd x ]nZ^ 2 \x,-x\y,-y) . (4.61) 

y >x 

Taking y — > x on the rhs of (4.59) with Nf = Nj, = 1 thus leads to 

xd x xGf\x) = 2x(z$(x)z£ 1 \x)+zW(x)z(+ 1 \x)) , (4.62) 

in complete consistency with (3.37) derived from the replica limit of the Toda lattice equation. This 
immediate consistency comes about since the supersymmetric generating functionals satisfy exactly 
the same Toda lattice equation as the fermionic and bosonic hierarchy of partition functions. 



5 Consistency Conditions and Toda Lattice Equations 

It is interesting to note that the Toda equations for QCD4 are special cases of a more general relation 
which was derived in ref. [14], and there referred to as "Consistency Condition II" due to this equation's 
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origin as a consistency condition in Random Matrix Theory [21]. This more general relation reads, 
after rotating to real masses, 

^ /+2 WW>) = \ Nf) (5-63) 

(x - y ) z v (0*/}) 

TV/ 

Z^ JV/+1) (y,{/x / })(^ nfdp, + x5 a; )^ iV/+1) (x,{^ / }) - (a; y) 

Taking the completely mass-degenerate case x = y = [if for all /, and switching to the conventional 
normalization for the mass-degenerate case, we recover the Toda lattice eq. (2.17). Similarly, if we 
take x = y = fif for / = 1, . . . , Nf and [if = z for / = Nf + 1, . . . , Nf + n we recover the Toda eq. 
(2.18). It appears that the single equation (5.63) is the generator of all Toda lattice equations. 

In ref. [21] it was stated that an analogous consistency condition exists for the QCD3 case, but no 
details were given. Here we supply this missing equation, which turns out to be of a quite different 
kind. We begin by noting that the kernel K(x, y, {/i/}) associated with the Random Matrix Theory 
(3.48) can be written in terms of the even-flavored effective QCD3 partition function [21] 



N f 



- zf Nf+2 \ix,iy,{m}) 



/ 4 2JV,) ({/*}) 

while in general from Random Matrix Theory it is known that it also can be represented in terms of 
the associated orthogonal polynomials P n (x,{[ii}) (see e.g., ref. [19] for details), 



f v x y 

(5.65) 

For this particular Random Matrix Theory ensemble the orthogonal polynomials split into two disjoint 
sectors of odd and even order. Now, since also these orthogonal polynomials in the large- N limit can 
be expressed directly in terms of the finite- volume partition functions [21], 



P2n(,x,{m}) 



^ Nl+1) (ix,{^i}) 



P 2N+ i(x,{fii}) = — (5.66) 

we can combine eqs. (5.64)-(5.66) into a consistency condition that must be satisfied by the effective 
partition functions. In contrast to the QCD4 case [21] we do not need to expand to first non-trivial 
order in 1/N, and the relation will therefore be algebraic. We find, after rotating into real masses 
throughout and fixing the normalization constant, 



zf Nf+2 \x y {»■}) = 



x — ' 



The three main types of partition functions that enter QCD3 are thus not independent. 



(5.67) 
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A particularly useful special case of eq. (5.67) is that of y = —x. Using Z^ Ns+1 \— x, {hi}) = 

—Z^L Ni+1 \x,{fii}) as well as Z^ Nf+1 \— x, {[J>i}) = Z^ f+1 \x, {pi}) we get a simple factorization 
identity, 

^N f+ 2 )( . -, , 2 Z^ N ^\x,{^})zf N ^\x,{^}) 

(x,-x,{ni}) = --2 6 . (5.68) 



X 



This special case actually also follows from combining the relations of Theorems I-II, and re-expressing 
results in terms of QCD3 partition functions alone. However, the more general relation (5.67) seems 
to go beyond what is contained in these theorems. 

The relation (5.68) gives rise to non-trivial identities among spectral correlation functions if we combine 
it with the replica method in the same way we derived the general relations (3.55)-(3.57). We can 
illustrate this very easily by focusing on the spectral density itself. To this end, consider the special 
case where we add 2N fermions of paired (common up to a sign) masses ±fi to the 2Nf paired masses 
Pi. The resolvents of the different theories can then be related: 

Gf Nf+2 \x, -x, {m}, n) = hrn ^<9 M \nzf Nf+2N+2 \x, -x,{m},/j,) 



lim — T 

N^o 2N 



df, In Z ^ + Nf+2N+1) ( x , {fr}, n) + df, In z ^ N f +2N+1 ) ( X) {^j, //) 
-d,lnZf N ^ N \ { ^,)] 
= Gf^ +1 \x, + G (2Nf+1 \x, - G? N '\M,ri (5.69) 

which in turn implies an identity among the spectral densities, 

2 [4 2 ^ +2) (A,z,-z,{M) +/4 2iV/) (A,{M)] = ^ +1) (A,x,{ ft }) + ^ +1) (-A,x ) { W }) 
+p ( i Nf+1 \\,x,{^}) +P l N f +1 \-X, x ,{^}) . (5.70) 

Needless to say, this general relation also follows from combining eqs. (3.55)-(3.57). 



6 Conclusions 

The purpose of this paper was to demonstrate that the replica method based on Toda lattice equations 
successfully can be applied to QCD3 in the e-regime. The result is a series of exact statements about 
spectral correlation functions for the Dirac operator of that theory in the phases where spontaneous 
symmetry breaking occurs. The Toda lattice equations were derived on the basis of known Toda 
equations for the effective QCD4 partition functions in the e-regime and exact relations between the 
two theories. Our results give the detailed support for general relations that formally can be derived 
by the replica method applied directly on the relations between QCD3 and QCD4. Because our 
corresponding Random Matrix Theory ensemble incorporate a number of determinants, these results 
go much beyond, but include as a special case, what has already been shown for the ordinary unitary 
ensemble. 

Armed with a reliable tool for computing spectral correlation functions by means of the replica method, 
we have also attacked a case that previously has not been considered in detail in the literature: an 
odd number of flavors in a regularization with an odd number of Dirac operator eigenvalues. We have 
shown that the peculiarities surrounding this case do not lead to pathologies in the spectral properties 
of the theory. This can also be understood from the point of view of the corresponding Random 
Matrix Theory once the singularity at vanishing quark mass has been treated carefully. 
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As in 4-dimensional QCD, the replica method based on Toda lattice equations allows for a smooth 
connection to the supersymmetric approach. In particular we have conjectured the form of the super- 
symmetric generalization of the relationship between the effective QCD3 and QCD4 partition func- 
tions in the extreme limit of the e-regime. With the corresponding 4-dimensional expressions known in 
closed analytical form this gives us a very compact expression for the analogous expression in QCD3, 
an expression which will be extremely tedious and difficult to derive directly in the chiral Lagrangian 
form (previous results in that direction do not include the so-called Efetov-Wegner "boundary" terms 
[23]). We suspect that a simple proof may be obtained by reverting to the Random Matrix Theory 
expression for the partition function. 

Finally we have illustrated how Toda lattice equations for QCD4 follow as special cases of a more 
general expression derived originally from the Random Matrix Theory representation. Interestingly, 
when we extend the proof in the most straightforward manner to QCD3 we do not recover our Toda 
lattice equations, but instead an algebraic relationship between the involved effective partition func- 
tions. Again, as a special case this relation can be reduced to an equation that can also be derived 
on the basis of the connection between QCD3 and QCD4. When applying the replica method to this 
relation, now justified on account of the Toda approach, we obtain exact relations between the spectral 
correlation functions involved. We have verified explicitly in simple cases that this general constraint 
indeed is satisfied. 
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